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In the former paper on this subject I have dealt with the formation of Laplace's 
dynamical equation for the tides, and the integration of it, subject to the limitation 
that the solutions obtained should be symmetrical with respect to the axis of 
rotation. In the present paper I propose to extend the method of solution so as to 
free it from this restriction. 

The diflficulties experienced by Laplace in his attempts to integrate the equation 
in question were so great that he abandoned all efforts to obtain a general solution, 
and confined his discussion to a few of the special cases which present the greatest 
interest from a practical point of view ; even in these simple cases however his 
original attempts to express the solutions by means of the coefficients associated 
with his nam e were discarded in favour of series proceeding according to powers of 
a certain variable used to define the position of a point on the Earth's surface. 
These power-series have been further employed by Lord Kelvin''^ to obtain a more 
general solution of the problem, but the results obtained, though of considerable 
analytical interest, do not lend themselves well to a numerical discussion. Both 
Amyt and Kelvin condemn the employment of the surface-harmonic functions 
as inappropriate, but a profound conviction that the efforts of Laplace, though 
unsuccessful, were well directed, has led me to take up the problem again from his 
point of view ; with what success will be seen hereafter. 

I was originally led to attack the problem by a totally different method from that 
of Laplace based on the work of PoincabeJ and Bryan §, and the principal 
analytical results, both in this paper and in the preceding, were at first obtained by 

* *' On the General Integration of Lapiace's Differential Equation of tlie Tides.'' * Phil. Mag.,' 1875. 
t ' Encjc. Metropolitana.' Art. " Tides and Waves," Section III., § 116. 

% " Sur I'eqnilibre d'une masse fluide animee d'un mouvement de rotation." * Acta Math.,' vol. 7, 
p. 355, et seq. 

§ "The Waves on a Rotating Liquid Spheroid of Finite Ellipticity." * Phil. Trans./ A, 1889. 
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a very lengthy analysis similar to that used by the latter writer. The comparative 
simplicity of these results seemed, however^ to point to the fact that they might be 
more easily obtained by less pretentious means. The deduction of the formulae in 
the former paper from the differential equation of Laplace presented no serious 
difficulties, but in attempting to apply a like method to obtain the more general 
formulae of the present paper, I found that formidable obstacles had to be overcome. 
The method of integration now adopted seems to leave little to be desired for sim- 
plicity, considering its generality, but the fact that it has been built up partly by 
working forwards from the differential equation, and partly by working backwards 
from the results, must account for the apparent artificiality of the procedure. 

In all previous attempts at the solution of the dynamical equations for the tides, 
the integration has been effected by assuming that the expression for the tide-height 
could be expressed by an infinite series of terms of known form associated with 
undetermined numerical coefficients. The differential equations then lead to a 
difference-relation between a certain number of these coefficients from which their 
numerical values are to be evaluated. The numerical determination of the coefficients 
will be facilitated when this difference-relation contains as few terms as possible. 
Now it is found in the present paper that, without imposing any restriction on the 
period of the disturbing force, if the form we assign to the terms of the series for the 
tide-height is that of the tesseral harmonics or Laplace's functions, a linear relation 
involving three successive coefficients only may be deduced, provided that the law of 
depth is such that both the internal and external surfaces of the ocean are spheroids 
of revolution about the polar axis. This however appears to be the most general law 
of depth which can be employed without obtaining more than three successive 
coefficients in the linear relation in question, and consequently our discussion deals 
only with cases where the law of depth is subject to this limitation. 

In § 1 I have collected the principal properties of the functions used in the 
analysis. These properties are for the most part well known, but in consideration of 
the want of agreement in the notation employed by different writers, I have thought 
it best to briefly prove such of them as are required in preference to giving references 
to places where they may be found. Moreover I have thus been enabled to write 
the results in the exact form required for subsequent application. 

§1 2-4 deal with the integration of the differential equations and the deduction of 
the linear equations (31), (40) connecting the coefficients in the expansion of the 
tide-height. These equations, the analogy of which with the equations (23), (23a) 
of Part I. will be at once apparent, constitute the chief analytical results of the - 
paper, and the remainder is occupied with the application of these formulae to 
the discussion of the free and forced vibrations on lines similar to those adopted 
in Part L 

§§ 5-11 treat of the free oscillations^ the discussion being confined to the case 
where the depth is uniform. A period-equation is obtained, and an approximate 
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method of determining the higher roots is given. The approximations will not 
however be sufficiently close for the earlier roots, and consequently it is necessary 
to evaluate these earlier roots by trial and error. The method of procedure is 
indicated by numerical examples, and several of the more important i*oots are 
tabulated for four different depths of the ocean. The most interesting result is the 
existence of a second class of free oscillations besides those whose existence may be 
at once inferred by analogy from the simpler problem of the oscillations of an ocean 
covering a non-rotating globe. The characteristics of the oscillations of this class 
are discussed in § 11. 

In § 12 a general analytical solution of the problem of the forced vibrations due to 
any disturbing force is given, but as the analytical expressions obtained are too 
intricate to afford much indication of the nature of the forced tides, the various 
types of oscillation which occur on the earth are afterwards treated numerically. 

In certain cases, intimately associated with those actually occurring, the analytical 
expressions however admit of considerable reductions. These cases are discussed 
in § 14, where theorems due to Laplace and Professor Darwin are obtained and 
generalized. 

§§ 15-18 contain numerical examples of the evaluation of the semi-diurnal and 
diurnal tidal constituents. The arithmetic is considerably simplified when the 
period of the disturbing force is rigorously equal to half a sidereal day or a sidereal 
day, and consequently these cases are first dealt with and the results compared 
with those of Laplace. Additional examples are however also given to illustrate 
the effects of the departure of the periods from exact coincidence with half a sidereal 
day and a sidereal day respectively, the cases selected for investigation correspond ing 
with the leading lunar constituents. 

§ 1. Properties of Tesseral Harmonics. 

Let P;, {ii) denote the zonal harmonic of order n, Then P,, is the solution which 
remains finite when ft = ± 1 of the differential equation 



Let 



dl 



3 



(1 - li?) -^'1 + w (n + 1) P„ = (1). 

V'M = {i-t.^f^% (2). 



dfj; 

Then, on differentiating the equation (l) s times, we obtain 

(i - r) ^^TTi - 2 (^ + 1) /^ -^^ + (n - s) (ft + 5 + 1) ^ = 0, 
or, 

(1 _ ^J,.^) |i ^( I _ ^.y,. p. ^ „ 2 (. + 1) /i- (1 - l^'r ^ {(1 - ^r'' P'n] 

+ {n-s){n + s+ 1) P: = 0, 
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which on reduction gives 



d 
dfjL 



< 



(1 - ^') t 



+ i ^'^ {u + l) 



s 



1 



P 



!^' 







. (3). 



Thus P^ is a sokition of the equation (3) ; the form (2) shows that it does not 
become infinite when /x = ^ 1^ while from (3) we see that the two functions 
P', cos s^, Pjj, sin 6'^, or what is equivalent^ the two functions V^^e^'^^, are spherical 
surface-harmonics of order 7i, In our subsequent work the latter forms involving the 
imaginary exponential will be more convenient than the real trigonometrical forms. 
We shall therefore describe the functions P', (/x)e^''* as the tesseral harmonics of 
order n and rank 6\ In some cases it may be convenient to apply the same nomen- 
clature to the /^associated function" P,^ (/x), but whenever this is done, it must be 
understood that an exponential factor is implied^ though not expressed* 

The tesseral harmonics of course include as special cases the zonal harmonics 
obtained by putting 5=0, and. the sectorial harmonics obtained by putting s = n, 
while, in accordance with the definition (2)^ for values of s greater than ri we may 
suppose that ^l{ii) = 0. 

The principal properties of the tesseral harmonics which we shall require may be 
derived from those of the zonal harmonics. Thus, if we differentiate s times the well- 
known relation 

(?i + 1 ) P,,+i — (2n + 1) jaP,, + nP,,_i = , . . , / . (4), 
we obtain 






(2n + l) /x 






#'^iR 






d'V 



«-l 



LvUj 



0, 



which, on making use of the formula 



c/P 



M+1 



dl' 



gives 



dfi 

d^F 
^ • ^ dpf 



dfjb 






df? 



df? 



(2r6 + 1)/^ '-~ + ('/I "4" s) 



«-l 



0. 



(5), 



On multiplying by the factor (1 ~ ju,~)»* this may be written 

(« - s + 1) P;.-i ~ (2n + 1) iiPl -{- (n + s) P:;„i = 



« » 



(6). 



Again by diflerentiating the equation (2) we find 



dV^ 



5/, (1 _ i^^)^^^ + (1 -t^'"-^ ^~~^ 



dfjb' 



SflPi + (1 -^ /X^) 



ms 



cP 



dfjj 

d 



s+l 



dfjif \y ^ ' dfjL 



stj:Pi -^ (1 »^ ij?f 






n (?^ + 1 ) 



d'-^T 



n 



duf ^ 



s (s 



\ tt^ X fi 



n 



25U -y~~ 4-5(5=-™ 1 ) - 
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in virtue of the differential equation (1) for P„ ; and therefore by naeans of (5) 






- ,„P _ (^^ - .9 + 1) (^ + g) ,p, _ p, . 

which with the aid of (6) may be expressed in the form 



Let us write for brevity 



^^ /^x .# 1 



D = (1 -— /x'^) , 



I ' t • • • • \ /• 



'^=---<»)-ii.==r--„i (!*'-■'') 



Then the equation (3) may be written 



1 

w***' 



while, if (T denote any constant quantity^ we obtain from (6), (7) 

(D + .,) P. = - '^^'f^ p;„ + '^-i-^;^^ifii> P;,-. . (10). 

The relation between the operators D, A may be written in the forms 

(D - cr/.) (D + cr/.) - (s^ - crV) .= (1 _ ^2) (A + cr) 1 

(D + <r/.)(D - (T/i) - (5^ ~ oV) = (1 - /*') (A - 0-) I ■ ' ^^^' 

which will be useful hereafter, 

§ 2. Transformation of the Dynamical Equations for the Tides, 

The formation of the differential equations for the tidal oscillations of the ocean 
has been fully dealt with in Part 1. It is there shown {§ 4) that, if U, V denote the 
northward and eastward velocity-components in latitude sin""ifi and longitude ^, 
when the system is executing a simple harmonic vibration in period 27r/X, these 
velocity-components will be expressible in terms of a single function t/; by means of 
the equations 



a (X3 _ 4^y ) a^ a y" (1 - fi^) (X^ - 4o) V) d<l> 
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Supposing that U, Y, \p are each proportional to e'^^^"^^^\ we rnaj take dxp/dcl) r=z isxjj^ 
and therefore, if we put 2a)s/X •=^ cr and introduce the abridged notation of the 
previous section, we may write the above equations m the form — - 



V'(l-^^)U = 

y(i-F)v = 



%scr 



1 



2(0 a s^ — a^/jL^ 



(J) — afi) xfj 



■^ 



2o)a 



1 



s^ — cr^/ji^ 



fiT) 



S'^ 



— fx (D — cr/x) 



O-^/L' 



1 

(T 



xfj ' ^ 



r 



(13). 



'/' 



J 



Tlie equation of continuity is 



3r 



1 



■|v(i.-f.V'U} + | 



;?.¥ 




(1 - /^^) 



where h denotes the depth, and Z the height of the surface-waves. On substituting 
for U, V from (13) and performing the difierentiations with regard to ?^, ; ^ this 
becomes "^ 



(i-/^^)c 



4co^a^ 



D 



h 



.9^ 



a^/uL^ 



(D —- cr/x)i// 



■ BM^i -l rt 



o"- 



4oj%^^ 



, V /T-) 



<T /X' 



(T/x) \p 



4(o^a 



2^^ 



or 



^'(]--/x^)^=(D + a-/.) 



A 



,1 



o-^/x"^ 



-^ ^J) _ 0-^^ ^ .«« jixjj . 



(14), 



This equation is equivalent to the equation (17) of Part T. A second equation for 
the determination of the two functions xjj, { is obtained from the pressure-condition 
at the free surface. On reference to § 2 of Part I., this condition is seen to lead to 



^ 



V 



gC + y 



, . (15), 



where v denotes the surface-value of the potential due to the harmonic inequalities, 
and V the surface-value of the disturbing potential. 

In order to effect the integration of these equations, we introduce two auxiliary 
functions ^i, '^g? connected with t// by the relation 



x}j = (D + a-ii) ^1 + {s^ - o-V^) ^ 



2 



(16). 



On applying the operator (D — cr/x) to the two members of this equation, we 
obtain in virtue of (11) 

(D - o-/x) r// = (1 - ix^) (A + 0-) ^i + {s^ - aY) % 

+ (^3 _ ^2^2) (J) _ ^^) qr^ _ 2crV (1 - /x') ^3 . . (17). 



Now the functions ^^, "^^ have as yet been subjected only to the single condition 
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(16). We may therefore impose on them any other arbitrary condition not 
inconsistent with the former. Suppose we choose them so as to satisfy the relation 

{A + (r)^^ = 2(T^lJL% (18), 

The two equations (16), (18) serve for the complete definition of the two functions 
"i^i, ^.2 ? making use of the latter, (17) reduces to 



(D «- cTfi) xjj = {s^ - aY) {^1 + (D 



cr/^) %} 



• ft 



(19). 



Thus on replacing i^ by its value in terms of "i^i^ ^^ in the right-hand member 
of (14), we deduce 



4S> <J 



( I - /x^) C = (D + a-fi) {h^, + /^ (D - cijl) %] ^h [(D + ail) H^^ + {s^ - crV ) ^ J 



If we suppose that h is constant, the terjns involving ^i will disappear, while in 
virtue of (11) we shall obtain 



^"^-f{i-i^')t, = h{i-i^'){^ 



a' 



O")^ 



3 



or 



(^--)^^ = ^^ 



(20). 



We have now for the determination of the functions i/;, C, ^\, ^^ the four 
simultaneous differential equations (15), (16), (18), (20). 



3. Integration in Series of Tessercil Harmonics, 



Let us suppose that i//, 4 v, ^F^, ^^ are each expressible as series of tesseral 
harmonics of the same rank s. Omitting the exponential factor e^'^'^'^^ we assume 
that 



7i •— OO 



^l>= t nP;(/.), 



n = s 

11 z- CO 



b — s ^,X % (/^)> 



n = s 
■n = CO 



v= t yf,P; (ix), }■ . 



• • . e ♦ e 



n — s 

n — CO 



^¥,= % <p:,(/.), 



n ~ s 
n = 00 



%= X ^K (/-), 



11 — s 



{ jLj 1. j t 



Then, if p denote the density of the water, and ctq the mean density of the whole 

\ Oli. CXCI. — A. XJ 
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system inclusive of the ocean, by well-known properties of surface-harmonics it 
may be shown that 



^Pff 



{M. +l)crQ 

Thus from (15), on replacing the quantities involved by means of their expansions 
in terms of associated functions, we obtain 



^1 ifct^ii V/^/ 



gm. 1 



3p 



(271+ l)(r(- 



P^ (/.) + Srf.P; iv) ; 



whence, if we equate coefficients of P'l in the two members, we deduce 



n=~ 9>A^l + yi 



909 09«98 



where we have written for brevity 



9u = 9 i 1 



o 

op 



(2n + 1) ctq 



9 e e * c 



* « 



From (16), (18), we have 

i// = (D + cr^) ^Fi + 6'2>?3 - ^ix (A 



4" cr) ^ 



S S » « 9 



But by means of (10) we find, on replacing ^F^ by its expansion. 



(D + (T/x) ^1 =:: —» S 



ot 



II 



\n — g-) {n — cs + 1 ) p, 

2n + 1 






while from (9), (6), we obtain 

/x (A 4- cr) ^\ ~ S< (cr — ii (ti 4- 1)} 






j^ 



n + 1 



uBu- 



s 



Thus 
(D + cr/x) 



I/. (A + a-)] ^Pl 

{/l(?Z' — 1) + O-} (76 

2?'t -f 1 



2^^n 



S 



+ 1) 



(22) 



(23). 



(24). 






and tlie right-hand member of (24) is therefore equal to 



t^^^] 



5 



+ iS< 



(,, (,, ^ 1) -f e^j (,, ^ g + 1) {(r^ +1) (y. + 2)_+ aUii +j) 

271 + r J-.-H-1-+ 2 76-M "■^' 



Hence, on comparing the coefficients of P', io the two members of (24), we obtain 

r^ ^ ,m^ , 1 (n--s){(n^l){n^2) + a} , (n + s + l) {(n+2)(n±^)±a} 



27Z,- 1 



276 -f O 
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Again, on replacing ^j, "^^ by their expansions in (18), the two members may be 
expressed as series of associated functions by means of (9), (6). Thus we find 



- S< [n (n + 1) - a-} P^, = 2(7^2^:, 






n ■+■ s 
2n + 1 






whence, on equating coefficients, we deduce that 



< [n{n+ 1) — a] 



2(y' 



01 — s 



2n 



I Pn~-\ "T i 



71+5+1 



2?i + 3 



Hn + l 



• a 



(26). 



Finally from (20), on expressing the two members by their expansions and 
equating coefficients, we obtain 

(n(w+l) + o-};S?, = -^Ct (27). 

The relations (26), (27), enable us to eliminate the auxiliary constants ^J„ aj,.i, 
<^.i from (25), and thus to express Tl as a linear function of Gj^^o? Q> Qi+2' On 
substituting for <^i, <^.i from the formula (26) in (25), we obtain 



n = 5^A - o- 



3 0^ - s) {(71 -- l)(n-2) + a} 
(2^- l){(n^ I)??.-- a] 



''-'-' ^i-. + ^.:^ 13^. 



2n - 3 



cr 



3 



(n + g + 1) {(n + 2) 0^ + 3) + a} 
(2n + 3) Jin + 1) (?2,~T2)"^^^^^~oy" 



^ •— s + 1 



2?z. + 



r ^^'' + 1 



2?^ + 1 



n + s + 2 



2n + 5 



A +2 



, 0^ -^ g) (^ - ,g -^ 1) {(^ ^ 1) (,, - 2) + a} 

(j .^ ^. .^ ^. , . ^x , Pn—2 



(2n - 1) (2n - 3) {(n - 1) n - o-} 

I r 2 _ 3 0t- 8) (^+s){0t-l)(^-2) + o-} _ g (w-g-f l)(?^ + g + l){(7^+2)(^>^ 4 3)4-q^}' 
"•"L^ "" (2ti-l)(2%+l){(«,-l>i-o-} °" {2ii+l){27i+3){{n+l)(n + 2)-crf 



|8:. 



O" 



(^^ + g + 1) {n + s + 2) {(n + 2) (?^ + 3) + (7} 
{2n + 3) (2?i + 5) {(71 + 1) (^ + 2) - a} 



^ 



n+2i 



and this, by means of (27), gives 



hV 



n 



{n — s) (n •— s — ' 1) 



4coV (2n - 1) (27^ - 3) {(n •- l)n - cr} 



Of^_2 "" A^^GJ 



+ 77 



{n -\- s ■■{- 1) ('?i -4- s + 2) 



(2^^ + 3) {2n + 5) {{n + 1) (7^ + 2) -- cr} 



p.9 



9 • 



where 



A' = 



S^^ 



{n — s){n -{- s) {{n — 1) {n - 2) + a] 



(28): 



'■ii 



cr^ {?i (?z + 1) + a) (2qi — 1) (2n + 1) {(7? — l)7^ — cr} {n (n + 1) + o-} 

(n + s -i- 1) (?z. — s + 1) {(71 + 2) (n + 3) + a-} 



(2?i + 1) {2n + 3) {{n + 1) (71 + 2) — cr} {n (n + 1) + cr}' 

V 3 
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This expression for A' may be somewhat simplified if we separate it into its com 
ponent partial fractions ; we thus find 



A! 



s'^ (n — s) (n -f s) 

<T^ {n{n + 1) -f^} """ {2n - 1) (2?i -f 1) 



n^ {n (n — 1) ~ cr} n^{7i (n + 1) + cr} 






(ti — s + 1) (n 4- s 4" 1) 
~~~(2n+ 1) {271 + 3) 



(^ + 2)3 



KK^**** 



2ti + 3 



(i?i + 1)3 {(71+ l)(?i + 2) — cr} (?2-fl)3{n (?i4-1) + <t} 



1 



{?Z. (?1 + I) -f- 0"} 



... j ,. 



(^^ 4- 1)3 _ s3 



_^ 0-3 ' # (2ti 4- I) {n 4- 1)3 (2?i 4- 1)„ 



(71 — 1)3 (?Z. — s) (71 4- <^) 



(?i 4- 2)3 (7^ — s 4" 1) (^ 4- s + 1) 



7i3 (2'/^ - 1) (2^^ 4- 1 ) {n (^i - 1) -- a} (/i 4- 1)^ (2?i + 1) {2n 4- 3) {(?i 4- 1) (n 4- 2) - cr} ' 
whence finally, remembering that cr = 2ms/X^ 



it 



X3 71 (7^ 4- 1) — 2ft)s/X 

4^ Ti^^^Tiy^ 



(^/i — 1)3 (n — s) (# 4- s) 



7^3 (2^1 "-• 1) (2??. 4- 1) {{n - 1) ?i ~ 2ft>s/X} 

(71 4" 2)3 (72/— s 4- 1) jn 4-^4-1) 
(^rTI)^^M^ - 2mslX} 



s \ iw «7 I • 



The relation (28) will hold for all values of n equal to or greater than s^ provided 
we suppose that CL2 = and C;.-i = 0. If we put for brevity 



iAjr, 



n 



% 



(71 -- s 4- 1) (^ — s + 2) 
{2n 4- 1) (2?i 4" 3) {0^ 4- 1)0^^) - 2^Js/X}' 

(fi- 4- s 4" 1) 0^ 4-^4-2) 
(2*1 4- 3) (271 4- 5) {(71 4- 1) (^^^ 4- 2) - 2fi)s/X} 



> . . . (30)j 



it may be written 



JiF. 



n 



4ffi3^2 



— ^n-^^S-Ai-t "^ -^iiMi + 3/B^*rt+2* 



R63placing T' by its value in terms of C:^, y^ [equation (22)]^ we obtain 



av^_2^tt-2 "~ ^n^n "T Vn^i 



+ 2 



4ct)3a3 



a e c 9 ;- 



where x;'„25 ?/» ai'^ defined by (30) and 



J-j 



^n 



4a)3a 



-, + A::, . 



e o 



e 9 



(31). 



(32), 



A', being defined by the equation (29). 

On putting 5 = it may readily be verified that the equation (31) reduces to the 
equation (23) of Part I. The manner in which such an equation may be utilized for 
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the determination of the free and forced vibrations has been fully discussed in that 
paper* A strictly analogous procedure might be adopted in the present case, but it is 
found convenient to modify the previous treatment in some respects. We shall 
therefore only indicate briefly the course of procedure when our analysis corresponds 
with that given in Part I., giving greater detail as regards those points where a 
different method has been found desirable. 



§ 4. JExtension to the Case of Variable Depth, 

The formulae developed in the preceding sections apply only to the case where the 
depth is uniform. We may however obtain a relation of the same nature as the 
equation (31) when the depth is a function of the latitude given by the formula 



h 



k + l{l'-fi^) . (33), 



that is, when both the internal and external surfaces of the ocean are spheroids of 
i^e volution. 

The above expression for h may be written in the form 



h = K -^ I 



4ft)2 



/^ 



2 



where 



K = k + l(l 



K + 



X3 



l(s^ — a^/jL^) 



C7^ 



4co^ 



(34), 



(35). 



Substituting the expression (34) for h in the equation (14) we obtain 



o,.o. 



4:(o^a 



(T 



3 



(l-/.^)C=(D + o-/x) 



K 



.9- — a'^fi 



i72 (^ — O-^) '/' 



K\JJ 



+ ^ [(D + cr/.) (D - cr/.) - {s^ - aY)] rp, 



which, with the aid of (11), becomes 



40)2^3 



cr 



(1 ^ ^^) 



^ 



I 



4:(jd^a^ 



(A — cr) i// 



(D + o-j^) 



fC 



(T 



a^fj? 



(D — (Tfx) Xp 



— KXJJ, 



The right-hand member is of the same form as that of the equation (14), except 
that h is replaced by k. We may therefore introduce two auxiliary functions ^i, 
%, defined by (16), (18), and proceed as in § 2, and we shall obtain in place of (20) 
the equation 



{^-o-)% = 



4:C0hc^ 



a^/c 



C- 






(A-cr)^ 



• • • » I ou )• 
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But J if we introduce the expansions (21) for xjj.^ 4 we find 



4.co'^xi' 



^ " lii (^ - '^) ^ = ^ CJ, + T-^ {n {n + l) + <T] Tl 



^mhi? 



P' 

X. r,t 



(37), 



whence the final equation which replaces (31) will be obtained by replacing C' by 



<^« + rx;i (^^ (^^ ■+ i) + ^}r«, 



4(»%i 



« » » 



(38) 



in the right-hand member of (28), and A by k in the left-hand member. 
Thus we obtain 



k\ 



Is 

n 



4<^(^ 



/Y'^ 

^n-2 



I 



CU + -, ^ {{n - 2){n -^ 1) + cr} Tl 



Am^a^ 



A' 



+ ?y« 



I 



C'?^ + ^r^" (^ (^' + 1) + o-} rf, 



I 



C-+2 + y:j:, {{n + 2) (n + 3) + o-} ^.^ 



4:m^a 



« 9 



(39) ; 



and therefore^ on separating out the parts of VI, due to Cjj y% respecti%^elyj we find 



where 



and 






« « e 



!.•§ M i l i- M /^^^ 

b'4 """^ '^it 



1 



^^« 



4:0)%^ 



n{^% '\' 1) -^ 



2ms' 



Vn 



Vn 



i-r3^{0^ + 2)(^i + 3) + '"^ 



4©%s 



X 



^j^ '°™~* 



4 » 



"4" A' 



s 

n 



i-ife^'^^'^ + i^+T 



r}.« _ 



l l— | w i w 



I 



4ft)^a^ 

K 



(n ~ 2) (?2 — 1) 4" 



2 £505 



^^i^-By^-^a 



e 9 



>■ 



*^ 






4 I 



in iwnn 



^ 



4{»"a^ 



'cos 



n (n + 1) + -^ !> Ai 



$ 



,8 « ,.S 



('/J + 2) (n + 3) + -v" !-2/;r«+2 



2a>s 



(40), 



e ? 9 



(41)= 



. (42). 



For the special case s = the equation (40) reduces to the equation (23a) of 
A. art) x» 



§ 5. The Period-equation for the Free Oscillations in cmi Ocean of Uniform Depth, 
To determine the periods of free oscillation we may proceed exactly as in § 6 of 
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Part I., making use of the equation (31) instead of the simpler equation (23) of 
the previous paper. On putting y", zero, the equation (31) gives 

an equation which must hold for all values of 7i equal to or greater than s^ it being 
understood that Q^2 = 0, and Cj_i == 0. The series of equations typified by (43) 
may be divided into two groups, in the former of which the suffixes involved are 
such that n-- s is even and in the latter odd. The types of motion resulting from 
these two groups may be treated independently, the former being characterised by 
symmetry with respect to the equator and the latter by asymmetry. The treatment 
of the two groups of equations will be exactly similar, and we shall therefore in the 
main confine our discussion to the former group. 
If we introduce the notation 



H 



K 



^^ 


^»^W 








^^l "^ 


71 ~ 


LI - 


— 




• » * 


u 






_2 


T s 


-• . . . ctd 




il 






inf. .. 



<" • • . ♦ « « (4:ttK 



it may be shown as in Part I. that provided U Qj = 0, 



'/^=oo 



f^g — ~ -*•■*-/{, 3 ^^^, — -^^,1 • • • • « » t y^O Jp 

and therefore the equation (43) may be written 

whence the period-equation for the free oscillations of symmetrical type is obtainable 
in the form 

when ^^ ■— 5 is an even integer. 

The same equation will apply to the asymmetrical types if we suppose that n — 5 is 
an odd integer, and that the continued fraction H^ terminates with the partial 

quotient '%f^ , 

In particular, putting n — s = 0, we can express the period-equation for the 
symmetrical types in the form 



'''siJs '^54 2^8+2 



^^ [/ — 1/ -^ acl v)if """ • . . . * (46a), 

while, putting n — ^* = 1, that for the asymmetrical types may be written 

'"**^ L^+s — 1^5+5 — . . . Ctd inf, ^ \ /• 

On the analogy of the problem dealt with in Part I., we may anticipate that, when 
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X lias a value in the neighbourhood of a root of the equation Lj, = 0, the continued 

fractions H'^„9, K'.+o will rapidly converge to small values. Further5 with large 

values of n^ the numerical values of Hj,„2> K',+2 tend to become equal with opposite 

signs. Hence there will be roots of the equation (46) which approximate to roots of 

the equation 

T/ — (47^ 

Let us therefore examine the nature of the roots of this latter equation ; this 
may best be done by considering the graph of the function AJ,. Putting y = A',, 
X = X/o>^ Ave have to consider the form of the curve 



2/ = 



•^ ^^ n (n 4- 1) — 25/^* 



4 



X 



n^(n + 1)^ 



(n - ly (n - 6') (7h + s) _^^ 

n^ {2n - 1) (2^T^l)Jn(n~'^l) -^ 2s/4 

(n -f 2y(n — s -f 1) (n -f s + 1) 



(n + 1)' (2n + 1) (2'?i + 3) {(n + 1) (n + 2) - 2six} 



It is evident that this curve will have two rectilinear asymptotes parallel to the 

axis of 2/, whose equations are 

2s 25 



X 



n (n — 1) 



X 



(n + 1) (n + 2) 



and that on passing these critical values, with increasing x^ the sign of y will change 
from positive to negative. The curve will pass through the origin, while when x is 
very large it will approximate to the parabola 

^n(n 4- 1) — 2s /x 



y 



4 *^ 



n^ (n + 1)^ 



Hence it must consist of three branches as in the annexed diagram^ where the 
dotted lines represent the rectilinear and parabolic asymptotes. 




,.2r 
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The roots of the equation Lj, = 0, regarded as an equation for the determination 
of X/oj, will be the absciss83 of the points of intersection of this curve with the line 

Since h is essentially positive/ the roots will all be real, and they will lie in the 

intervals between 

2s 2s 

+ 00 . 



— 00 



(n -f 1) (71+ 2)' n{n-l)' 



For large values of lig^jA^iJ^o? the two extreme roots will approximate to the roots 
of the equation 



while the remaining two roots will approximate to 



2s 



2s 



{n -h l)(n-h 2)' n(n-l)' 

The two former roots are those which have their analogue in the special case 
treated of in Part I., for which 5 = 0, in which case they have equal magnitudes but 
opposite signs. These roots, we may expect, will approximate to roots of the period- 
equation, at least when n is large. 

In order to see the signijficance of the remaining roots, it is convenient to transform 
the period-equation into a diflferent form, which moreover is far better adapted for 
the more accurate numei'ical determination of the earlier roots. 

§ 6. Modified Form of the Period -Equation. 



Referring back to the equations (29), (30), (32), v/hich define <, yl, L', we see 
that (31) may be written in the form 



(n—s) 



{2n-l) {n(n'-l) -2(os/\} 
\2 n(n-i-l) — 2ft)s/X 









Of^ 



+ 



4«2 n^ (n + If 

(n+s+l) r ^n+2Y(n-s+l ) 0^4-5+2) ' 

(2n+3){{n + l)(n + 2)-2o)s/X}\ (n + l)\2n + l) '' "^ (2^ + 5) "+^ 



hj 



g^s 



4a)%: 



(48). 



Hence, if we introduce a new set of auxiliary constants, D^, Dj+i, &c., such that 
for values of n equal to or greater than s 



(n 4-1)^ {n — 
VOL. CXCI. — A. 



^ns , ^' (n + s + 1) p^ 



n {n + 1) 



2cos 



D 



n 



(49). 



X 
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the equation (48) reduces to 



2?i 



(n + If ( ^~s) j^, n^n-h s-i- 1) ^^ 

" X^ f , , ^ . 2u)S 



n^ (ft + 1) 



3 _% 



?i 



4ft>^a^ 



C'l + '/^^ (^ + 1)^ : 



Hi 



4a)^a^ 



(50) 



Thus, if we write for brevity 



M; 



4^ 



n {n + 1) 



. 2a)S 



n^ (n +- 1) 



2 



hgn 



4ce)%^ 



N;:^?^,(7^ + l) 



2(08 

X 



>■ 



» « » 



-^ 



and put 7,t = 0, the equation (43) is replaced by the two following : 



(51), 



(n f If (n ^ s) 

2n ^ 1 '-^^ 



(n + 1)^ (n - g) ^^ 
2?^ — 1 



?i-~i 



N-' D' _L '"'^ (^ + g + 1) p, 



51 + 1 







= 



^ « • {OjLi)» 



For the determination of the symmetrical types we therefore have the series of 
equations 



M; Gi + - 



s^(2s + 1) ,^, 



LlS ~J- i 



3 



D 



s+1 "~" ^? 



0, 



(s + 2)M (g+ 1)^ 2^ + 2 ) _ 



{s + 3)1 2 , , (s + 2)H2s_+ 3) ,, _ 



e e 



On eliminating the quantities C% Dj+i, C'+g/ &c.v by means of a continued 
fraction we find the period-equation in the form 



MJ 



a 



a 



s + l 



a: 



'S±2 Q 

N"::4-i - M::+a - N:^.^ ~ . . . ad inf. "^ 



« 5 



(53), 



s+l 



where we have written for brevity a\^ in place of 

n^ (n + 2)^ (71 — s + 1) (t^ 4- s + 1) 



a fi ft 



© e 



• (54). 



In like manner the period-equation for the asymmetrical types may be written 



N 



al 



a 



s+l 



W 



s+2 



M 



s-fl 



isjs _ ]v/r 



s+3 



, , , ad inf. 







6 f » ♦ 9 



(55). 
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We may also write these equations in a variety of alternative forms in which 
prominence is given to any one we please of the quantities M*, Nj,. These forms 
are obtained by giving n different integral values in the equations 



M 



N 



n 



s 
91 



_^^ n~l 


^■M-2 


- n;_s - 

- Mf,_8 - 


• t 




„ w+l 


- MJ,,, - 


- M?,^3 - 


• • • 




■^^ n-~2 


^.^1 - 


• * • 



= 0. (56), 



. (57). 



In each case the former continued fraction terminates with a partial quotient 
involving a\ in the numerator and either Mj or N,' in the denominator, while the 
latter proceeds to infinity. 

For the symmetrical types, if we use the form (56) we must suppose n — s an 
even integer, whereas if we employ (57) n — s must be supposed odd. The reverse 
will of course be the case for the asymmetrical types. 

The continued fractions of the present section will not converge so rapidly as those 
of the preceding, but in spite of this drawback they present considerable advan- 
tages. In the first place the numerators of the partial quotients, which are obtained 
by giving n different integral values in the expression (54), are independent of X. 
These numerators, which further are in a convenient form for logarithmic computation, 
may therefore be tabulated once for all, v^hereas the numerators of the partial 
quotients in the continued fractions of the last section require to be re-determined at 
each successive trial in attempting to solve the period-equation by trial and error. 
Moreover, the evaluation of the denominators Mj„ N;\ by means of the forrnula3 (51) 
may be very quickly effected, even though a fairly large number of these denomi- 
nators is required, whereas the evaluation of the quantities L^ by means of (29) and 
(32) is extremely laborious. 

Another disadvantage resulting from the use of the preceding form is that when 
X/co is near the value 2sln{n +1) the functions xl_i, yf,_i both become large, while 
L-,_i, L',^.1 have both a zero and an infinity in the immediate proximity of this 
value. Hence, in order to evaluate Lj,«i, Lj,+i for a value of X in this region it is 
necessary to observe a very high degree of accuracy in the numerical work. The 
singularities which occur in the left-hand member of (46) when X passes through one 
of these critical values no longer appear if we write the period-equation in the form 
(57) with the value of n appropriately chosen. 



§ 7. Expressions for the Velocity-Components. 

The auxiliary constants Dj, Dj+j, &c., introduced in the last section, may be made 
use of to express the velocity-components by means of series of surfiice-harmonics. 

X 2 
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Thus from the first of equations (13) we have^ on replacing \jj by its value in 
terms of "iTj^ "¥^^ 

^ (1 - /x^) U = - ^^ [^Fl + (D ~ cr^) ^4 



On introducing the expansions for '^j, ^Fg on the right we find 



.v/(l~/x3)U 



isa 

2o)a 



a: 



n 






4-l)(?H~6') 



2n + 1 



135 



tVl I al 



n ^^n 



(n i- a — I) (n — s) (7 1 ~ <r + 2) (11 + ^ + 1) ^, 

Pit— 1 I o,, , o P/i + l 



2?i - 1 



2^^ + 3 



which, with the help of (26)5 (27), reduces to 



.y/(l ■— JUt^)U 

or from (49)5 



2misct 



a 



h 



2P 



(?^ — .9)(^i + q^ — 1) p, , (w. + s + 1) Oi ~ (t) p, 



n+ 1 






?l 



2ti — • 1 






mI 



2n + 5 



J- 91 • e ^ O ) • 



This may also be written in the form 



/ 



iXa 



^ (i_^3)u = ^yS{/.ci--D:„}p;. 



• s ^ « 



Again, from (13 



v/(l-/.2)V = 



'l(T 



IJL^(1 — JUt^) U 4" 



(T 

2wa 



xb 



h 



aX 



. (59). 



+ ^ fit ifiCi ^ Bt:) Fi + :^ tr^j^i \ . . (60), 



whence, by means of (6), we may express ^{l -- pf^)Y by a series of surface- 
harmonics. 

The corresponding formute when the depth is variable may be obtained by 
replacing h by k and Cj, by 



Ci + j:jz^^ \ n {n + 1) + -^ I ■^'^* 



4ft)"a' 



X 



-'- W3 



so that we find 



V^(l -"/x^)U 



'iXa 



V 



fC 



jtl 



<^'^+ii«^^^^^+l)+^?}^' 



"15 

n 



-^ 



Dl IP 



>j, / -*- p^ 



% 



v/(l-)a^)Y = ^^Sti- 



/C 



^^+da^{''^'+^'>+'? 



K61) 









J 



where the quantities D', are now defined by 
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Qi-1 4* 



+ 



n^ (n 4- s + 1) 



,,(,,^i) + 2^^jr:.„, 



2?i + 3 



I r / . -, \ / . ^v . 2cl)s 

— — — j' 



C:.,. + ^^-, j(n + J) 0^ + 2) + "YJI^Ux 



(62). 



The formulae (61) cease to be of use in a special case which will present itself here- 
after for which k = 0, It will be seen in a later section that the expressions on the 
right become indeterminate in this case, so • that the determination of the velocity- 
components must be effected by means of the formulae (12) or (13). These latter 
formulae seem at first, sight to indicate that the velocity-components become infinite 
in latitude sin~^ (c^/cr), but the forms (61) indicate that such cannot be the case, at 
least when k is different from zero. 

§ 8. Approximate Determination of the Higher Roots of the Period-Equation, 

If we take the period-equation in the form (56), and as a first approximation omit 
the continued fractions from the left-hand member, it reduces to the quadratic 

Ml = 0, 
or 

^^ f / . .\ 2«sl ., , .„ hg. 



4(» 



2 



{»(«+ 1) -^} -■«'('> + i)'i5ii=o. 



For large values of n the roots of this equation will give a sufficiently accurate 
approximation to the roots of the period-equation, since it may be seen that the 
continued fractions tend to limits comparable with |-#, and therefore small in com- 

parison with n^ {n -f 1)^7~^ , when n is very large and \/a> has as its value either of 

the roots of this equation. We may even obtain a fair approximation by omitting 
the term containing 5, in which case the formula for X corresponds with that obtained 
when the rotation is omitted. 

A better approximation will however be obtained by representing the continued 
fractions by their first convergents instead of entirely neglecting them. The 
approximate form of the period-equation is then 



or 



4ft)3 n^n + If #(2n - 1) (2^i + 1) {{n -• l)n - 2o}s/X} ' 

{n -f 2)^ {n — s -{ - l){n -f- ^ + 1) ^ hgn _ 

In -f If {2% + l)(2'?^ -i- 3) {{nVl)in + 2) - ^wsjX} '^ 4. ^ ^ "~ 






We thus get back to the equation I4 = 0. The roots of this equation may be 
approximated to numerically by Hoener'b process, the significant roots being those 
which lie in the intervals between •— cx) and 0, and between 2sjn (n — 1) and + ^ • 
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For the particular case s ^= 0, the biquadratic to which the equation Ll = is 
equivalent reduces to a quadratic, the two roots which remain finite being of equal 
magnitude and opposite sign. This special case has been examined in Part I., and it 
will be seen on reference to the tables there given (§§ 7-8), that the roots of the 
equation L^, = give a very good approximation to the roots of the period-equation 
except in the case of the earlier roots when hgjior'a^ is smalL In the present paper 
I have examined in some detail the special cases corresponding to the values 1 and 2 
for s^ and the approximation is found to be equally rapid, as will be seen from the 
tables given hereafter. Consequently, in these cases at least, all except the two or 
three smallest roots will be obtained with adequate accuracy by finding the roots 
which lie in the stated intervals of the equations L', = with difl'erent integral 
values of n. 

The roots so found will not however form the complete series of roots of the 
period-equation. We may in fact anticipate that the remaining roots of the equation 
L^^ z=z will also approximate to roots of the period-equation. To obtain a better 
approximation of the roots of this class, it will however be preferable to make use 
of the period-equation in the form (57). As a first approximation we omit the 
continued fractions and obtain 

Nl = or n {n + I) — "^^ =i 0. 

This method of approximation will be valid if when X/w = 2s/ii{n'^ l) the two 
continued fractions involved in (57) are small in comparison with n {n + !)• B^t it 
may readily be verified that with large values of n these continued fractions become 
comparable with co^a^/hg, and therefore the desired condition will certainly be 
satisfied when n is sufficiently large. 

A better approximation may be obtained by representing the continued fractions 
by their first convergents. We thus obtain as the approximate form of the period- 

equation for the determination of the root which lies near 7 , :. 



or 



(n — 1)2 (n + If (n — s) (n + s) 



, , ^ . '2ws (2?^ — 1) (2^ -I- 1) 



y*y<ffTW'yifp!" 't^ 



4«M ^ ^ X 






71^ ill + 2)3 {n ~ s + 1) ill + s 4- 1) 






2^; 



{n + 1) {n + 2) 



(2% + 1) (2^1 + 3) 

aMIMII I UMlU ' JW M MMI I I MI MWM a iUW ' J I " ■ 

2 ms ] ('?i + 1 J' ((ti -f 2)2 lign^-\ 






X 4a)%2 



ANALYSIS TO THE DYNAMICAL THEORY Of THE TIDES. 159 

Using the first approximation in the terms on the right, we deduce 

(n - lY(n -f 1)^ (n - s) (n -\- s) n^n + 2y (n -g+ l)(n + s + 1) 

, , ^ . 2a)S ~ (2n - 1) (2n + 1) ^ (2n + 1) (2n + 3) . ' 

w(ti+ 1 ]"---—=== — — — — — — — y— — — (63). 



This formula is found to lead to the roots of the period-equation with a surprising 
degree of accuracy. 

Our analysis is only applicable when h is small in comparison with a, but subject to 
this limitation the approximations of the present section will improve as hg/Aa)\i^ 
increases, that is, as the depth of the water increases or the angular velocity of rotation 
diminishes. They will give good results even with small values of 7h when c«> is 
sufficiently small, and they may be used to determine the limiting values assumed 
by the roots when the angular velocity of rotation is indefinitely reduced. 

We see then that^ the roots of the period-equation are of two classes, which may 
be distinguished by their limiting forms when the rotation is annulled. The roots of 
the former class are such that the values of X remain finite when cj = 0, their limiting 
values being given by the formula 



± /s/ 



ii (n + l)hg 

cf 



There will be an equal number of positive and negative roots of this class, but 
though these approach the same limiting values their numerical values will not be 
equal as in the case where 5 = 0, and the positive and negative roots must therefore 
be determined independently. 

The roots of the second class are all positive and are such that the values of \/o) 
tend to finite limits when o) is reduced to zero, the limiting values being given by 
the formula 

(o n(n + 1) 

whereas X will tend to the limit zero. 

The analogue of the types of motion which correspond with the former roots will 
still be oscillatory when the rotation is annulled, but the types of motion corre- 
sponding with the roots of the second class will cease to exist as oscillations when 
the angular velocity of rotation is reduced to zero. These types of motion will 
have their equivalent in steady motions, but an infinitesimal amount of rotation 
would immediately convert such steady motions into oscillatory motions of very 
long period. 

For the particular case s = the roots of the second class are all zero even 
when the angular velocity of rotation is finite. Hence steady motions can exist 
on a rotating globe, but these are necessarily of zonal type. We have in fact 
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already seen in Part I/'" that the only forms of steady motion which, can exist are of 
this character, and have explained the fact by stating that the steady motions not 
of zonal type which can exist on a globe wdthout rotation most have their analogue 
in the more general case in oscillatory motions whose period bears a finite ratio to 
the rotation-period 5 no matter how great the latter may be. Our present wTjrk 
confirms this statement and throws further light on the nature of these oscillatory 
motions. 

§ 9. Evaluation of the Earlier Roots. 

The errors resulting from the use of the approximate formulae of the last section 
may be considerable in the case of the earlier roots for w^hich n has small values. 
To obtain these earlier roots w^e must therefore proceed by trial and error, the 
preceding method being made use of to obtain values with which to commence 
the trials. 

As a concrete example we will discuss in detail the computation of the positive 
root of the first class corresponding to the case n = 4^ s = 1, when the depth is 
given by hg/ia)^^ = yq. Taking p/o-g = 0*18093, and introducing the numerical 
values of n^ s, and A, the equation L| =: becomes 

{k/o)y^ - 0-3333 {\/ojf - 5-5481 {X/o)f + 1-0906 (X/o)) - 0*0418 ~ 0. 

By HoRNEB^s process the greatest positive root of this equation is found to be 

2-43265. 

Now experience shows that the numerical value thus suggested is in general too 

small.t We therefore select for a first trial a value rather larger than that 

indicated, say, for example, 

k/o) = 2-4400. 

From the formula (54) we find 

log al = 0-2553, log al — 1'1652, log al — 1*7289, 
log al = 21450, log a? = 2-4769, log al = 2-7537, 
log aj = 2-9915, log a^ = B '2001, log aj = 3-3859, 

while the values of the expression 

^v" {n + If ^^'^'' 



for the values 2, 4, 6, 8, 10 of ^ are 

1-6046, 18-794, 84*517, 250^92, 589*4. 

t Compare the 2nd and 3rd cokimns of Tables I, and II., Part I. 
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Thus we find from the formulae (51), with X/w = 2'4400, 

M^= 6-104, Ml =9-753, MJ =- 23-25, MJ =- 144-97, Ml 
NJ = 1-180, N^= 11-180, NJ = 29-180, NJ = 55-180, NJ 



- 426-9, 
89-18. 



It will be convenient for us now to introduce the following abridged notation :- 



ei = 









cv 



n-l 



a 



n 



a: 



W. + 1 



*n. 






M* — 






^w-1 


< 


^^?j + l 








< 




^n-2 




MU 


< 


<-l 


^71-2 



, , , ad inf. 



. . . ad inf. 






■^^ 7l~2 



* * .* 



(64), 



the last two continued fractions terminating with the partial quotient which involves 
al in the numerator, and either M^ or Nj in the denominator. 
From these definitions of the quantities e, f, E, F, we have : 



a 



^n — 






n~l 



^n J 71+1 



a 



71 






P == 






a 



^ 



n—l 



•*~^7l 



e 



71+1 



a 



71 



N' 



n 



^'7l~l 



(65) ; 



while the period- equation may be written in the forms : 

N;, - E:,_i - eUr = 



(66). 



Suppose that we neglect fu ', making use of the numerical values obtained above 
for the quantities M, N, a, by successive applications of the formulae (65) we obtain 



loge{o = wO-7556, log/g' = 1-2229, log e^ = nO-7829, 
log// = 0-9666, log ej = 710-9649, log/' = 0-5606. 

In like manner, if we neglect e\o, we find 

log/i = 1-2498, log e^ = 710-7801, log// = 0-9668, 
logej = nO-9649, log/' = 0-5606. 

Now the two values of /' obtained by these methods are respectively the 6th and 



yoL. cxci. — A. 



Y 
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5th convergents of the continued fraction f^. Since we find that to the degree of 
accuracy retained these are equal, it follows that all subsequent convergents are 
sensibly equal to either of them. Hence the infinite continued fraction f^ may be 
replaced by its fifth convergent without sensible error. 
Similarly we find 

log F} == 0-1834, log EJ — 0-5045, log Y\ — 0*8266, 
and therefore 

Ml — FJ - // = 9753 - 6707 - 3-636 — - 0*590. 

As a second trial we take 

X/cu = 2-4600. 

Proceeding as before, we deduce 

Ml- FJ-Zs'^: 10*234 - 6-619- 3-607 = 0-008. 

We conclude that there is a root of the period-equation lying between 2*4400 and 
2*4600 ; by interpolation its value is found to be 

2-4597. 

The same method may be used for the determination ot the roots of the second 
class, the initial trial values being suggested by the formula (63). As a numerical 
example, if we put ^i = 5, 5 = 1, hgjico^ar = yo ^"^ (^3)? we find 

2co/X=: 40-974, or X/co := 0'04881. 

For a first trial we take 2a)/X ::=^ 41, and deduce 

Ni — El - 4 = — 11 + 8-678 + 2-602 = 0*280. 

As a second trial we take 2o)/X = 41*280, and obtain 

W,-El- el- - 11*280 + 8-657 +2-593 = - 0*030, 

and therefore, by interpolation, 

NJ- El --6^ = 0, 
when 

2co/X =: 41*253, or X/co = 0*04848. 

I have selected for special investigation the asymmetrical types when s :=^1, and 
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the symmetrical types when 8=2, these types presenting special interest hi relation 
to the diurnal and semi-diurnal forced tides. The annexed tables s^ive the values of 
the computed roots, together with the corresponding periods of free oscillation for 
the four depths of 7260, 14,520, 29^040, 58,080 feet treated of in Part I., corre- 
sponding to the values -^q, ^o, yq, and ^ for %/4oi%^. 



L— Depth 7260 feet {hg/Amht^ = 4^, p/o-Q = 018093). 





( 


s = l. 

..-, J^ 






5 = 2. 

A 






'» 


' 




"''*"' ~"- s 


Approximate 
root. 


Corrected 
root. 


Period* 


Approximate 
root. 


Corrected 
root. 


Period. 


• 

hH 

3 

• 

hH 
hH 

W 

1— 1 



Wr=:2| 

w- = 4< 


1-8331 

- 1-0906 
2-0042 

- 1-8472 

0-6932 

0-05702 

0-03853 


1-6337 

- 0-9834 
2-0685 

- 1-8234 


hrs. mins. 
14 41 
24 24 
11 36 

13 10 

• 


1-4745 

- 06582 
2-0024 

- 1-7006 


1-3347 
- 0-6221 

1-9866 
--- 1-6595 


hrs. mins. 

17 59 
38 34 
12 5 

14 28 


n = l 

n = 3 

n = 5 


0-6401 

0-06251 

0-03784 


- 

days hrs, 

1 13 

16 

26 10 


0-1532 

0-08091 


0-1671 

0-08178 

■ 


days hrs. 

6 
12 5 



II.~Dbpth 14,520 feet {hgjioy^a^ = -^q, p/a-Q = 0-18093). 











_^,_ . . 




s = 2. 

, .. . . A 1 








"^ ^ "" "" ^ 


f 




' "■ "^ ^ 


Approximate 
root. 


Corrected 
root. 


Period. 


Approximate 
root. 


Corrected 
root. 


Period. 


* 










hrs. mins. 






hrs. ndns. 


hH 
m 


n = 


.i{ 


1'9588 
- 1-3021 


1-8677 
1-2450 


12 51 
19 16 


1-6804 
0-9128 


1-6133 

- 0-8922 


14 52 

26 54 


Q 


'it — 


1 


2-4327 


2-4597 


9 45 


2-4363 


2-4349 


9 51 


hH 
hi 

02 


?i — 


- 2-2949 


- 2-2907 


10 29 


- 2-1706 


2-1547 


11 8 


0-7502 


0-7283 


davs hrs. 

i 9 


• • 


• • 


days hrs. 

• • 




n = 


<> 

o 


0-08423 


0-08673 


11 13 


0-2054 


0-2129 


4 17 


Q 


n — 


5 


0-04881 


0-04848 


20 15 


0-10067 


0-10081 


.9 22 



=* Throughout these tables the periods are expressed in sidereal time. 
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III. Depth 29,040 feet 


(%/4a)%^ = 


^ /Ol P/O-Q - 


: 018093). 








s- 1. 




r 


' 6- - 2. 

_.A. 




r ■■■■ • " 






" ^ 


Approximate 
root. 


Corrected 
root. 


Period. 

lirs. mins. 
11 5 
14 50 

7 40 

8 1 

days hrs. 

1 5 

8 21 
17 19 


Approximate 
root, 

2-0241 

- 1-2960 
3-1295 

- 2-8911 


Corrected 
root. 


Period. 


M 

m 
m 

1 — 1 

o 


n — 4i< 


2-2027 
™ 1-6439 

3-1183 
- 2-9958 


2-1641 

16170 

3-1274 

- 2-9961 


1-9968 

- 1-2855 
3-1293 

— 2-8856 

0-2554 
0-11472 


hrs. mins. 
12 1 
18 40 

7 40 

8 19 


Class II. 


w = 1 
n = 3 
n— B 


0-8213 
0-1116 
0-05636 


0-8149 
0-1124 
0-05625 


0-2523 
0-11472 


days hrs. 

• « 
3 22 
8. 17 
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.V. Depth 58,080 feet 


, [hg/io)^a^ = 


- i p/^o — 


0-18093). 










s~- 1. 






f = 2. 










r"" 


A 






j\^^ ._ 






™^ 


r- '-- ""■■■" " ^j 


Approximate 


Corrected 


Period. 


Approximate 


Corrected 


Period. 


— - — — -- 






root. 


root. 




root. 


root. 


hrs. mins. 










hrs. mins. 








-{ 


2-6431 


2-6288 ^ 


9 8 


2-5636 


2-5535 


9 24 


m 


n — 


- 2-1726 


- 2-1611 


11 6 


1-8623 


- 1-8575 


12 55 


m 

o 


n = 


-{ 


4-1626 
■»- 4-0501 


4-1659 
- 4-0506 


5 46 
5 56 


4-1820 
- 3-9610 


4-1818 
- 3-9592 


5 44 
• 6 4 












days hrs. 






days hrs. 


h-\ 


n = 


: 1 


0-8886 


0-8873 


1 3 





• « 


• « 


m 
m 


n — 


z 3 


0-13354 


0-13375 


7 11 


0-2865 


0-2876 


3 11 


Q 


n — 


= 5 


0-06108 


0-06105 


16 9 


0-12332 


0-12332 


8 3 



The approximate roots here given have been evaluated by the method of | 7, 
except in the case of the roots of the second class for n ■==■ 1^ s-==- 1, where, instead 
of replacing the continued fraction involved in the period-equation by its first con- 
vergent, I have made use of the second convergent, so that the approximate form of 
the period-equation from which this root is determined is 



or 



N} 



MJ 



OA 



M» 



a 



or, what is equivalent, Lg = Q. 



N} 













0, 
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By a comparison of the approximate values given in these tables with the true 
values we see that for extreme cases here tabulated, the error involved in the 
approximation does not amount to more than about 3 per cent., even with a depth as 
small as 7260 feet. We have here a justification of the statements made in the 
last section as to the approximation to the higher roots. 



10. Oscillations of the First Class. Determination of the Type. 



We shall describe as oscillations of the first class those whose periods remain finite 
when the rotation-period is indefinitely prolonged, that is, those for which the roots 
of the period-equation are of the first class. The types of motion whose periods 
become infinitely long with the rotation-period will be called oscillations of the 
second class. 

The determination of the type involves the determination of the constants C", 
Dj^_i, Q+2J &c. (supposing for convenience that we are dealing with symmetrical 
types). For this purpose we may either make use of the formulae of § 6, or we may 
make use of the formulao of § 5 for the determination of the C's, after which the D's 
must be computed from equation (49). The latter method will be closely analogous 
to that used in § 9 of Part I., but the former is the more convenient when the 
numerical determination of the constants is required. 

If we make use of the notation (64), the following relations may be deduced from 
the equations (52) :~ 



C 



2n + 1 






{n + s) (n - If 

2n + 1 



e! 



C 



n 



2n 4- 1 



n } 



Gl-i {n + $){%- 1) 



fs 



J)' 



{n 



C 



»+i 



• s + 1) {n + 2f 

2n + 1 

(w - s + 1) (%^+~2)2 



E 



"1 



> • • (67). 



F 



s 

n 



But we have seen in the last section how the quantities e, f E, F may be deter- 
mined numerically. Hence the above formulaa allow us to compute the ratios of the 
constants C, D. One of these constants must be regarded as arbitrary, and the 
ratios of the others to it can then be computed. When the type under examination 
is that which corresponds to a root of the period-equation approximating to a root of 
M^ ==: 0, we select as the arbitrary constant of integration the quantity C'„ as the 
continued fractions c, f E, F required to determine the ratios of the remaining 
constants to this one will then be free from singularities. 

When these ratios have been determined we may substitute their values in the 
formulse 



t =:: n^e*^^^"^^^ 



• • • I 



c 



g 



ns ns 

^n~^ ps , p5 , ^^4-2 

pg ^ 71—2. I ^ n I r^s 



-^71+2 ~r 



„ _^_±i_ ps 






■is 
n 



n-\- 



4- 
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y{l-/.^)U 






n' 



jJii , . . -f- 



^/(1~/^')V 



2a)a 



/x G',/^""' -^ '^'^ 



fi 



s 



ns ^i(\t + S(j>) 

aX ^^"^ 



C! 









- P^^ 4- P'^ 4. -^-^i±i. Ps _L. I 



J 



« • e 
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D 



■?i4-l 



» 



Oi 



p^l+1 + • 



?? 



« « 9 
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I 



D^' 



!„ —^LiL ps I 

• ' ' x" fig '*- /I— I I 



J)i 



ps -^ n+l 



H 



ns 



. . -j -—f- g^_^ Pn^2 + ff'H Ph + ■'~7^"5'« + 2 P71 + 2 + • • • 



C 



72. 



C 



fl 



I 



4. 

I « • * 



and we shall obtain expressions for the height of the surface-waves and the velocity- 
components. 

I have not thought it worth while to compute any of these series in detail, as the 
general character of them may be inferred from the series computed in Part I. for 
the special case where s = 0. For large values of n^ and even for comparatively 
small values of n when ligjioi^a^ is large, the quantities C will rapidly diminish as 
we pass away in either direction from C',. Hence the most important term in the 
series for ^ will be that involving Pj,, and this term will in general sufficiently 
predominate to decide the number and approximate position of the nodal parallels of 
latitude. 

If we neglect Cj^i in comparison with Cj+i and suppose that a» is small in com- 
parison with X, the formula (49) gives 

r (r + 5 + 1) QUi ~ (r + 1) (2r + 3) Dl 

and therefore Dj will be of the same order of magnitude as Q,+i, Hence, when r is 
less than n, D' will be of the same order of magnitude of C'+i, and similarly it may 
be seen that when r is greater than n^ .Dj, will be of the same order of magnitude as 
Q.„i. Thus the predominant terms in the expressions for the velocity-components 
will be those involving C",, Dj^-^i, D^^+j. 

The exponential factors indicate that the type of motion involved will consist of 
waves propagated round the sphere with uniform angular velocity X/6' about the 
polar axis, and that there will be s crests or troughs on each parallel. Positive 
values of X will correspond with waves propagated in the opposite direction to the 
rotation, that is westwards, while negative values will correspond with easterly 
waves. The paths of tlie fluid particles will be ellipses with their axes directed 
along the meridians and parallels. 

§ 11. Oscillations of the Second Class. 



In dealing with the oscillations of the second class we proceed in the same manner 
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as before, retainmg D^ as the arbitrary constant of integration when the type 
under consideration is that whose period is approximately given by the formula 

n (n + 1) — 2a}s/k = 0. 

For we may anticipate that this quantity will predominate over the others, at least 
when the depth is sufficiently large or the angular velocity of rotation sufficiently 
small. The ratios of the remaining constants to D'l may then be computed from 
the formulae (67), and on substituting these ratios in the equations 
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we obtain expressions for the height of the surface-waves and the velocity-com- 
ponents. The values of X for the oscillations of this class being always positive, the 
direction of the wave-propagation will always be westwards. 

By way of numerical illustration 1 have computed the series for ^ corresponding to 
the case n = 3, s = 1, and to the values i^, -^o, -^g, 4 for hg/Ao/a^. In these four 
cases the series within the square brackets are found to be 



2-977FI 



1-4735PJ 



0-7296Pi 



0-36l7P^ 



0-1880PH-0-3753P^-0-0916PJ+ 0-01153Plo 

- 0-00093PJ2-}- 0-000052P1, - 0-000002P}e -f- . . . 

- 0-3260P1 4- 0-11690P^ - 0-01309P> + O'OOOSOPlo 



- 0-000032P}2 + 0-000001 P 



u 



0-2248Pi + 0*03159PJ - 0-00168PJ + O'OOOOSlPJo 

-0-OOOOOlPk-f . .. 
0-1277P1 + 0-00814Pi - 0-00021PJ -f- 'OOOOOSPlo - 



It will be seen that as the depth increases or w diminishes the coefficients all 
become smaller, while the convergence of the series improves. It may be inferred as 
in the last section that CJ/DJ + , will tend towards the limit zero when r is less than 
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n, while Cj+i/D^ will tend towards a finite limit. In like manner when r is greater 
than n, Q/Dj^i will tend towards zero while C'„i/Dj. will tend to a finite limit* 

Let us examine the limiting forms assumed by 4 U, V, when the rotation is 
annulled. On putting 7i {n + 1) — 2a)s/\ = the relations (51) give 



L^ o>^M: 






U' 'El =z r (r + 1) *— ^^ (^ ^~ 1) =z (r — n) {r -f n + 1)- 



Thus from (52) we obtain, if r < n, 



If ^^ 


+1 


and if r > n^ 
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O.J, 
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(r + 1)2 (r "™ s) 



Hence if we retain only the most significant terms and put a>D| = A^, we find 
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If therefore we suppose that o) reduces to zero while A', remains finite, ^ will 
reduce to zero, but the velocity-components will tend to finite limits given by 
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in virtue of (7). 

Hence the steady motions to which the oscillations of the second class reduce 
when the rotation is annulled, involve no deformation of the free surface. This of 
course may readily be verified by a direct method. For if % v, w denote the 
velocity-components referred to fixed rectangular axes, it may be seen that all the 
conditions of the problem will be satisfied by 
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dv dy dy dy dy dy 

"^ oz dy dec az cy ^ oo:^ ' 

where ^ is an arbitrary function of x^ y^ % independent of t. These solutions make 
ux + t;^ + 1^2^ = Oj and therefore involve no deformation of the surface. If we refer 
them to polar co-ordinates, they are equivalent to 

^ sin ^ ac|) ^ B^ ^ ' 

or 

y(i - /x^)u = aj, v(i - ,.^) V = - (1 »- /.^)^^ 



8/i, 



This solution becomes identical with that found above if 






For the special case ^ = 0, the values of X corresponding with the roots of the 
second class will be zero, and the corresponding types of motion steady^ even when o) 
is finite. The types of steady motion will however involve a deformation of the 
free surface in all cases where the angular velocity of rotation is different from zero. 
These cases have been fully discussed in Part I., § 14. 

§ 12. Forced Oscillations, Geneml Analytical Solution, 

The problem of the forced oscillations involves the determination of the quantities 
Q in terms of y%^ from the equation (31) or (40). Dealing first with the case of 
uniform depth, and supposing that the disturbing potential involves only a single 
term y?^PJ^ (/x) e*^^*"^'*^^, we have to solve the simultaneous equations 



^jsMs "~" l-'5+2L/f-{-2 "T 2/sM-2^s+4 — t), 



•^H-2Q,^o — L?,Q + y],QU^ = hYJio>ht^, 



with the condition that U Cj = 0. 
From these we deduce, as in § 5, that 
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(r < n + 1). 



'^O 






VOL. OXCL—A Z 



170 



Al'* K^ I 



;. HOUGH OH THE APPLICATION OF HAEMOHIO 



Tims the ec| nation which involves yl^ Becomes 
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We may now deduce Cj_23 Qt^4? • • - Qi-yi^ Ct. 
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by means of the formute 
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It should be noticed that the term mvolving Pj, need not here be the predomi- 
nating term of the series within the square brackets ; it will however be so when 
the value of X for the disturbing force is in the neighbourhood of those roots of 
the period-equation which approximate to roots of L^, = 0, But if X have as its 
value another root of the period-equation, say, for example, one which approxi- 
mates to a root of Lj,, = 0, the series within square brackets will differ only by a 
constant multiplier from the series in the expression for the tide-height for the 
corresponding type of free oscillation^ since the equations which determine the ratios 
of the C*s are evidently the same in both cases. Hence^ for values of X in the neigh- 
bourhood of this one^ the predominating term will be that which involves Pj,. 
Consequently^ when n and m differ widely, the numerical computation of these 
series will become laborious; for as we proceed away from the term containing 
P^, towards that involving Pj^, the terms will at first increase in magnitude, and 
the convergence of the series will not assert itself until the term depending on 
P^,, has been passed. These circumstances will not however occur in any of the 
cases of more pi-actical interest. 

Whatever be the nature of the disturbing potential it will be possible to expand 
its surface- value in a series of surface-harmonics. Thus the most general value 
of V which can occur may be expressed in the form 



§2=:^ 51=00 



'^ "C* '^ 






The deformation of the surface due to each term may be calculated independently 



ANALYSIS TO THE DYNAMICAL THEORY OF THE TIDES. 



171 



and the results superposed, so that the deformation at tune t resulting from this dis- 
turbing force will be given bj 



c 



A 



4(»%3 



5 = 00 n — co 

5 = ll~S 



•^y(A(+^) ^ g.g-KA^ + s*) 



"LTs 



«--3 






+2 






Kfl 






t^-V-j — Q 



^^ 



The corresponding formulae when the depth is variable may be obtained by replacing 

h, yl, 8:^, HJ,, Kj,, Lj„ <, ^^., by k, G"^, Af„ Ht, 3^1, ll, 1^^, '^1 respectively, where 
K, Gj, H', ^^, 7j* are defined by the equations (35), (41), (42), and 



m 



n 



s 
n 



^nV 



% 






«-2 



^n~2Vn~2 ^nV'u 



% 



H 



nH-2 



• • • 



while A:^ is obtained by writing Si, in place of yj,. in the right-hand member of (42), 



§ 13. Classification of Tides. 

In the last section Ave have reduced the problem of the evaluation of the forced 
tides due to any disturbing force to that of the development of the disturbing poten- 
tial as a series of surface-harmonics. This development for the case of the disturbance 
of the ocean due to the attraction of the sun and moon has been already dealt with, 
and reference may be made to Professor Darwin's article in the ' Encyclopgedia 
Britannica' for a full account of it. We give here a short summary of the principal 
results of which we propose to make use. 

The principal part of the disturbing potential will consist of spherical harmonic 
functions of the second order, and when expressed by means of zonal and tesseral 
harmonics the terms which occur will be of three types, characterized by the rank of 
the harmonic involved. 

For the first type 5=0, and the corresponding tides will be expressible as series of 
zonal harmonics. For these types the value of X will be small in comparison with 
that of <y, so that the period of the disturbance is long compared with a sidereal day. 
The terms will cease to be oscillatory when the orbital motion of the disturbing body 
is neglected. The tides generated by these parts of the disturbing potential have 
been already dealt with in Part 1. 

The terms of the second type are those for which .9 = 1 ; in certain of these terms 
the value of X will be equal to o>, so that the period is rigorously equal to a sidereal 
day, while in the rest the period will reduce to a sidereal day when the orbital motion 
of the disturbing body is neglected. If n denote the mean orbital motion of the 
luminary, the '' speeds" of the principal diurnal tides will be co and &> — 2n, We 
propose to neglect the sun's orbital motion, so that for each of the principal solar 

z ^ 
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diurnal constituents we shall suppose that X = cu rigorously. The same analysis will 
then apply to one of the lunar diurnal constituents, while in order to illustrate the 
effect of the departure of the period from exact coincidence with a sidereal day we 
shall evaluate independently the lunar diurnal constituent for which 

K{=^ 0) -- 2n) = 0-92700 c^. 

The principal part of the tidal oscillations will be due to the third part of the 
disturbing potential, which involves harmonics of rank 2, The period for the tides 
due to these terms will differ but slightly from half a sidereal day, and will reduce to 
half a day exactly when the orbital motion of the disturbing body is neglected. We 
shall therefore assume that X = 2co rigorously for the solar semi-diurnal tides, while 
we shall take the value 

X/2co = 1 "^ n/o} = 0*96350 

as typical of the lunar semi-diurnal constituents. The analysis applied to the solar 
constituents will be rigorously applicable to the sidereal luni-solar semi-diurnal tide 
usually denoted by the symbol Kg. 



§ 14. Special Cases, 

Instead of making use of the equation (31) as we have done in § 12, we may of 
course compute the forced tides by means of the equations (49)3 (50) of §6^ 
determining incidentally the constants D*,. Thus, if we suppose that all the y's 
are zero except yj^-i, we have the following equations for the determination of D] 

_ N^D'^ 4- ^^ "^ "^~ C> 1 = 

(^±3)_^_2 {s + 2y (2s + 3) _ 



.<? 



where the terms on the right are all zero, except in the second equation. Now it is 
evident that if Nj •= 0, or 



X — ^"^ 



s -f 1 
all these equations will be satisfied if 



JJl = is+lf{2s+l) 



4ca%^ 

and 

^^S + l -"^ ^-^s + 2 -^ ^3 + B 
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It; follows that if the disturbmg potential be of order s -^ I and rank s, and the 
period be |- (^ + 1) days, the tide will involve no rise and fall at the free surface, but 
will consist merely of horizontal currents. If we put 5=1 the requisite period will 
be rigorously equal to a sidereal day, and the circumstances will correspond with those 
we have assumed to characterize the solar diurnal tides. We therefore conclude 
that in an ocean of uniform depth the solar diurnal tides will involve no iise and fall. 
We shall however see hereafter that for certain of the lunar d^iurnal tides the 
difference between the period and a sidereal day may be sufficiently great to render 
the rise and fall of considerable importance, unless the depth is very small. 

We have seen in § 4 that the formulae applicable to the case of variable depth may 
be deduced from those applicable to the case of uniform depth by replacing h by k 



and C;\ by Cl + 



n (n + 1) 4" -—■ 

A- 



■ps 



But if Yn = 0, rj == — gnQi, aiAd therefore when X = 2o)/{s + 1) and all the ys 
ave zero except yj^^i, 02 .is? Cj+s? &c., will all be zero, while 

GUi + 4^ {{s+l){s + 2) + .. {s + 1) } rUi = 0; 
whence we obtain 






1-2 (s -f 1)^ "^^ 



Thus the forced tide will be similar in type to the disturbing potential which 

produces it, though it will be inverted unless I < or > j—-~-- 

(s+ l)V.+i 

This theorem will admit of application to the solar diurnal tides on putting 5 = 1, 
in which case it reduces to a theorem given by Laplace.^"' The critical value of /, 
for a system comparable with the earth, is considerably greater than such depths as 
occur on the earth ; hence, for depths comparable with that of the ocean, the diurnal 
tides will be inverted when the ocean is deeper at the equator than at the poles : 
they will however be direct when I is negative, so that the ocean is deeper at the 
poles than at the equator. 

It is evident that the equations typified by (3l) will all be satisfied with the Cs 
all zero if ^ = 0, since in this case the right-hand members will reduce to zero. In 
like manner the corresponding equations which apply to an ocean of variable depth 
will all be satisfied when «: = 0, if for all values of n 



" "^ 4fA2 






. F = 



* Lamb, * Hydrodynamics,' § 212. 
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This equation leads to 



CI 



is 
n 



4:0} 



% 



n (n + 1) -f — >- 



X 



-fro/" I 



2o)8' 



There exists then a certain law of depth, depending on the period^ for which the 
tide will always be similar in type to the disturbing potential which produces it. 
This law of depth is expressed by the formula 



Iv 



3 



A= I 



If we suppose that X = 2<y rigorously, it reduces to 

so that the depth will be a maximum at the equator, and will gradually decrease on 
passing away from the equator to zero at the poles."^^ 

For other values of X the formula for h will make the depth negative at some parts 
of the surface unless I is positive and X > 2o}, The latter condition does not occur 
with any of the leading tidal constituents,t but it would hold good in the case of the 
semi-diurnal tides due to a satellite whose motion in its orbit was retrograde. If 
however we neglect the mutual attraction of the waters, the theorem under dis- 
cussion may be supposed to apply to an ocean covering that part of the surface over 
which h is positive, the remaining parts of the surface being supposed to consist of 
continents. When ^ is positive, these continents must, for the lunar semi-diurnal 
tides, reduce to small circumpolar islands, while for the same tides when I is nega- 
tive they will cover the whole globe with the exception of two small seas surrounding 
the poles. 

For the diurnal-tides, the shores most coincide with pai^allels of latitude approxi- 
mately 30° north and south of the equator, while for the tides of long period the 
appropriate forms of sea will be bounded by two parallels nearly coincident with the 
equator. A change in the sign of I in all cases involves an interchange between the 
seas and the land. 

It should be noticed that the formuhe (12) make IT, V infinite at the points where 

/x r- ^ — , that is, at the shores. This indicates that the neglect of the squares of 

the velocities is not allowable in the neighbourhood of the coasts no matter how small 
the amplitude of vibration may be, and seems to point to the existence of ^* breakers" 
as an essential accompaniment of the tides. 

* Cf, Lamb, * Hydrodynamics,' § 213. 

t There will be smaU tidal constituents depending on liiglier powers of the moon's parallax for 
which X exceeds 2a>. Of, Daewik, '* Harmonic Analysis of Tidal Observations.'' ' Brit. Assoc. Eeport,^ 
1883 (Southport), § 3. 
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I 15. Solar Semi-diurnal Tides. 

In the last section we have considered some special cases in which the tide-height 
is expressible by a single term. In general it will however only be expressible by a 
series of terms. It may be shown, as in § 5 of Part I.^ that this series will be finite 
when the law of depth is such that 






0,'' 



where n is an integer, n — ^ being even or odd according as we are dealing with the 
symmetrical or asymmetrical types of rank s. The values of I determined from this 
equation will in general involve X, but for large values oin^ they will approximate to 
the same values of I as those required for the expression in finite terms of the long- 
period tides, since for such values 2o)s/\ will be small compared with n {n + !)• 

In other cases the expression for the tide-height will involve infinite series. We 
deal in the present section with the case where I is zero, so that the depth is uniform. 

The numerical computation of the semi-diurnal tides admits of special simplicity 
when X = 2a) exactly. Putting X = 2co, .9 = 2, in the formulae (30), (29) we obtain 

^'' ^ (n + 3) (2n + ij (2n + 3) ' ^'^ "~" ^M + 3) {2n + 'E) ' 

a2 _ fa >- 1) (n + 2) (n ~~ If (n + 2) (n ^ 2^ (n ■-' 1) 

'' "^ n^ {n -f 1)'^ 7^3 (n + 1) (2n - 1) (2?i + 1) n (n + If {2n + 1) (2% + 3) ' 



the last of which gives on reduction 

A^ = 



2 _ 2(n-l) (n + 2) 

■n 



9i('?i + 1) (2'?i - l)(2^i + 3)' 



This general formula for A| fails to hold when n = 2 ; for, in this case, 7i — 5 and 
t^(n -— 1) — 2(os/X are both zero, and therefore the second fraction involved in the 
expression for A| is indeterminate. To determine it^ limiting form we must first 
suppose the period slightly different from half a day, so that X is not rigorously equal 
to 2(0 ; the formula (29) then gives 



X^ 2.3-4a>/X 4M.5 



^^ ^ 4a?2 2^3^ 33.5.7{3.4-.4ft)/X} ' 



which, on putting X = 2co, reduces to 



A2 = -zr-; 



2^33 3^7.10 5.7; 

The formulae for xl, yl, A| are now in a convenient form for logarithmic computa- 
tion, and we may readily deduce the following numerical values 

* Of, Laplace, ' Mec. Oel.,' Part 1., Book IV., § 7. 
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J 



Taking hgl^o^^a 



A^ 


= 0-085714 


Al: 


- 0-0233766 


A|: 


= 0-011544 


Ar 


=i 0-006823 


-og 


£c|= 3-75696 


log 


xl = 3-63639 


og 


a;|= 3-45469 


log 


xl— 3-29451 


log 


2/1 = 2-67778 


log 1/1= 2-14569 


logi/i — 3-81531 


\ogyl— 3-57512 




= iV , and p/o-Q = 


T2 . 

JU2 ' 


= + 0-063428 


U-- 


= - 0-0001156 


U-- 


- - 0-012412 


Lia ■ 


= — 0-017379 



a2 
A 2 



= 0-00449 
^ 0'003179 
^ 0-002367 
^ 0*001830 



1 



og XlQ 



log a;?3 

log ^14 



3-1564 
3-0360 

4-9297 



log ylo 
log y% 

log 2/?i 



d'dbbO 

3-0993 



0-18093, the formula (32) leads to 



T 2 _ 

ijio - 


z ^ 0-019860 


T 2 _. 


: - 0-02128 


T2 _ 


: — 0-02216 


T 2 _ 


: — 0-0228 



Thus if we neglect Kjs, and make use of the formula 



we obtain in succession 





K:.=: 




-2^w~-2 






3Ssion 


Mi 






logKf6 = 


=: r^ 5-671 






log K^ = 


>i3-03947 


log K?4 ~ 


= n 5-9226 






log K| = 


n3-72835 


log K?2 = 


=: 77.4-2165 






log K| - 


2-71588 


log Kio •= 


= //?. 4*5753 











and therefore, since C'+a/C* = K^^^^/yl, we find 



log (q/ci) 
log (c^/q) 

log (Cl/Cl) 
log(CVCi) 



0-03810 
nT-58266 
wi-22417 
?iT-0002 



log (cycfo) 

log (CyCf,) 

log (cycL) 



« 2-8300 
n2-6914 
n 2-572 
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But if ®f,P', (/x)6'^^^'^'*^ denote the height of the ^equilibrium' tide resulting from 
a disturbing potential y', P:J,, (/x) e^<^^-^''i'\ we have 



yl 



gMl 



Therefore 



c:; 






hgjicohc 



T2 Tr2 

JLJi> ^~' JL\, i 



rm 



and on introducing the numerical values for Kg^, LS, Kl we obtain 



Cl=- 1-94:76 m- 



Thus we have : 



log (Cy®i) = » 4-9646 
log (Cy®i) = 5-6560 
og (Cy©|) = « 6-228 ; 



log {Cy<El) = n 0-32760 
log{C^/@;f) = 1-91025 
hg{Cl/(Bl) = «T-13442 
log(Cy©l)= 2-1346 



and therefore if we suppose the exponential or trigonometrical factor to be involved 
in (El, so that the height of the equilibrium tide is expressed by QtlVl (fi), the height 
of the corresponding dynamical tide is given by 



^ = ©^ [- 1-9476 P:^ ~ 2-12617 P^ + 0-81331 Fl - 0-13628 P| 

+ 0-01363 P?o -~ 0-000922 Pf^ + 0-000045 P?, - 0-000002 Ff, + . . .' 



For points lying on the equator we have /* = 0, and it may be shown that in this 

case 

. ., 3.5 .. . (2?i 4- 1) 



J^2rt ( ) 



whence we deduce 



2.4 .. . (2'ft - 2) 



loo" Fi 



log F] 



loa- P« 



0-47712 
w 0-87506 
1-11810 



log P?o = 1-4325 



log Pj2 = n 1-5464 



WP1=: n 1-29419 



log Pi.i 



1-643 



logP?6 = nV73. 



Thus the values of the successive terms of the series within the square brackets 
for points at the equator are 

- 5-8428 + 15-9463 + 10-6746 + 2-6829 + 0*3690 + 0*0324 + 0-0020 + 0-0001, 
which, on addition, give 23*8645, But the height of the corresponding equilibrium- 
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tide at tlie equator is SQtl, and therefore the ratio of the height of the tide to that 
of the correspondiDg equilibrium tide at the equator is 

+ 7-9548. 

Ill like manner the tide-height in any other latitude may be compared with the 
equilibrium tide-height, but the process will be laborious in the absence of tables of 
the functions Pj,."^' 

The above example has been treated in some detail as ilhistrative of the method to 
be employed for the computation of the forced tides by infinite series ; the chief part 
of the labour is involved in the determination of the quantities oif^,, ?/' , A',, but when 
once these have been determined, since they do not involve the depth, it is easy 
without much additional labour to multiply cases for different depths. Besides the 
case already considered, which corresponds to a depth of about 7260 feet, 1 have 
computed the series for depths of 14,520, 29,040, and 58,080 feet, corresponding 
with the values i^q^ x^-, and -J for hg/iora'^. For these depths the series within the 
square brackets is replaced by 

- 0*83227P2' + 0-21694P;i - 0'02615P^ + 0'00180Ps' - 0-000080P?o 

+ 0*000003 P'i,+ . VM 

- 191'925P^+ 15'696Pj- 0-8082P5+ 0'0256P8' - O'OOOSPfo + . .. 
and 

1'9610P:^ - 0'06823Pa + 0'00164P;' - 0*000025?^ + . , . 

respectively, giving for the ratio of the tide-height to the equilibrium tide-height at 

the equator the values 

- 1*5016, — 234-87, + 2^1389. 

When the depth of the ocean is greater than 58,080 feet the tides are therefore 
direct at the equator. They gradually increase in magnitude as the depth decreases, 
and become infinite and change sign for some critical value of the depth rather in 
excess of 29,040 feet after which, for further decrease of the depth, they remain inverted 
until a second critical value is reached which is somewhat greater than 7260 feet, 
when a second change of sign occurs. The very large coefficients which appear when 
hg/ioy^a^ = Xq indicate that for this depth there is a period of free oscillation of semi- 
diurnal type whose period diff'ers but slightly from half-a-day. On reference to the 
tables of § 9 it will be seen that we have, in fact, evaluated this period as 12 hours 
1 minute, while for the case hg/4:0)^a^ = -^(^ we have found a period of 12 hours 
5 minutes. We see then that though, when the period of the forced oscillation differs 
from that of one of the types of free oscillation by as little as one minute, the forced 
tide may be nearly 250 times as great as the corresponding equilibrium tide, a 

^ The zonal harmonics from Pj up to P; Lave been tabulated by Glaishee, 'Brit. Assoc. Reports,' 
1879, but I do not know of tlie existence of any Tables of the Tesseral Harmonics, with the exception 
of a few giycn by Thomsok and Tait, 'Nat. Phil.,' vol. 2, § 784. 
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difference of 5 minutes between these periods will be sufficient to reduce the tide to 
less than ten times the corresponding equilibrium tide. It seems then that the tides 
will not tend to become excessively large unless there is very close agreement with 
the period of one of the free oscillations. 

The critical depths for which the forced tides here treated of become infinite are 
those for which a period of free oscillation coincides exactly with 12 hotu^s. They 
may be ascertained by putting X = 2a> in the period-equation for the free oscillations 
and treating this equation as an equation for the determination of A. The roots may 
be found by trial and error as in § 9^ the approximate values with which to commence 
the trials being suggested by the discussion already given. The two largest roots are 
found to be given by 

hg/ioy^a^ = 0'10049, lig/iaM = 0-02545, 

and the corresponding critical depths are about 29,182 feet and 7375 feet. 

We have hitherto supposed that p/ctq = 0*18093, but for purposes of comparison I 
have also examined the case where p/cr^ = 0, that is where the mutual attraction of 
the waters is neglected. The series for ^ in this case become 



k2 
.10 



1= eJi[l-0927PH"" l-91817Pi- 0-66909PS + 0-1070lP^ - 0-01036F 

+ 0-000683P?o - 0-000033P?4 + O'OOOOOIP^g - • • .] 

^ziz @2[-_ i.0733P',+ 0'24502Pi-.0-02790PH-0-00185P;--0-000080P?o 

+ 0*000002 Pro - . . .] 

i - m [9-34370P1 - 0-70311 Pa + 0-03449Pp^ - O-OOlOSP'g + 0*000022Pio - . . .] 
^ == Q^l [l-7739P<? - 0-05750P^ + 0-00132Ps - 0-00020P^ + . • •] 

or the depths 7260, 14,520, 29,040, 58,080 feet respectively. From these series we 
deduce as the ratio of the tide to the equilibrium tide at the equator the four values 

— 7-4343, - 1-8208, + 11*2595, + 1-9236, 

results which agree, except in the third case, with the numbers given by Professor 
Lamb^ deduced from the numerical formulae of Laplace. 

It will be seen that in three cases out of the four here considered the effect of the 
mutual gravitation of the waters is to increase the ratio of the tide to the equilibrium 
tide. In two of the cases the sign is also reversed. This of course results from 
the fact that, whereas when p/(Tq = 0-18093 one of the periods of free oscillation is 

* By a careful re-computation of the semi-diurnal tide for the case /3=10 (notation of Professor Lamb) 
I find tlie following series more accurate than that given for ^/H'" : — 

*^3 -f. 6-1915 1^^ + S'2U7p^ + 0-7234^8 ^ 0-0919 i^o 4- 0*0076 1^^3 + 0*0004 ».^4 ^ . _ 

This series reduces to 11*2595 when p =^1, thus agreeing with the result obtained above. 

2 A 2 
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rnther greater than 12 hoiirSj when p/cr^ = the corresponding period will be less 
than 12 liour 



I'D" 



§ 16« Lunar Semi-dmrnal Tides. 

A siniilai^ method to that of the last section may be used to evaluate the lunar 
semi-diurnal tides for which we take X/2c<j = 0*96350« The arithmetical w^ork is^ 
however, more severe, in consequence of the fact that the quantities r4, yl., A' must 
be evaluated from the formulae (29)^ (30) which do not assume the simple forms 
obtained in the last section. Substitutino^ in these fornuilEe the value of \l2o} 
quoted above we deduce 

A! = 0-075603 At,i = 0-003845 



Al = 0-01.9864 Af, =: 0-002750 

M — 0-009856 ATj — 0-002028 

Ai„ = 0-001968 ; 

log as'io = 3-1566 
log a4 = 3-0362 , 
log x{^ — 4-9299 ; 



Ai' = 


0-005834 


log X" 


3-76020 


log x\ — 


3-63756 


log a| = 


3-45530 


log Xg _ 


3-29488 


log 2/5 - 


2-68108 


log yi — 


2-14687 


log 2/1 = 


3-81592 


log yl = 


3-57549 



log 2/lo = 3-3867 
log y% = 3-2313 
log yii = 3"-0994. 



Our procedure is now exactly similar to that of the last section. Thus, if the 
height of the equilibrium tide be 

m Pi M, 

we find, when hg/Aoi^a^ = -^^ and /j/o-q = 0-18093, 

I - @| [0-10395 P| + 0-57998 Fj - 0-19273 P^ + 0*03054 P^ 

- 0-002960 P?„ + 0-000196 Pj, - 0-000010 P?4 +.. .]• 



Similarly, when hg/Ao)'^a^ = -^ 



•2 0; 



^ = ei [- 1-0647 P? + 0-24038 PI - 0-02774 P;-; + 0-001867 P's 

- 0-00U082 Pio + 0-000003 Pf, - 



a a » 



when hg/Aoy'xi^ — -- 



10> 



^ = ©I [9-1 181 P? -0-71533 PH- 0-03621 Pi 

— 0-001136 P^ + 0-000024 Pfo— . . .] ; 
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and when hg/iojp'a!^ ::=: \, 

C=:m [17646 P! - 0-06057 P^ + 0-001447 P^ - 0'000022 P^ + . . .]. 

From these series we find for the ratio of the tide-heights to the equilibrium tide- 
heights at the equator the four values 

- 2-4187, - 1-8000, + 11-0725, +1-9225. 

On comparison of these numbers with those obtained for the solar tides in the 
preceding section, we see that for a depth of 7260 feet the solar tides will be direct 
while the lunar tides will be inverted, the opposite being the case when the depth is 
29,040 feet. This is, of course, due to the fact that in each of these cases there is a 
period of free oscillation intermediate between twelve solar (or, more strictly, sidereal) 
hours and twelve lunar hours. The critical depths for which the lunar tides become 
infinite are found to be 26,044 feet and 6448 feet. 

Consequently this phenomenon will oocur if the depth of the ocean be between 
29,182 feet and 26,044 feet, or between 7375 feet and 6448 feet. An important 
consequence would be that for depths lying between these limits the usual pheno- 
mena of spring and neap tides would be reversed, the higher tides occurring when 
the moon is in quadrature, and the lower at new and full moon. ^^ 

There appears then to be a considerable range of depth comparable with the mean 
depth of the ocean over which the reversal of the spring and neap tide phenomena 
would take place^ but in that the actual tides are highest in the neighbourhood of 
new and full moon we conclude that the effective depth of the ocean does not lie 
within this range, and that none of the periods of free oscillation of the actual ocean 
lie between twelve solar hours and twelve lunar hours. The true effective depth is 
almost certainly less than 26,044 feet, and therefore both solar and lunar tides will 
be in the main inverted, though the configuration of the land and of the ocean bed 
will probably give rise to considerable variations of phase in different places. 

The shortest period of free oscillation of the second class for the case s=2 approxi- 
mates to, but is in excess of, three days. But if n denotes the moon's mean orbital 
motion, the speed of the lunar semi-diurnal tide is 

2 (o) — n). 

If we equate this to ^ o), we obtain 



n = f 0). 



Hence, if the moon's orbital motion were accelerated, or the earth's rotation 
retarded, until the month and day were in a ratio less than 6 : 5, it would be possible 
for the period of the lunar semi-diurnal tide to confound itself with one of the periods 
of the oscillations of the second class, and the tides would then tend to become very 
large. 

* Cf. Kelvin, 'Popular Lectures.' vol. 2, p. 22 (footnote). 
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I 17* Solar Se-mi-clmrnaLITides in an Ocea/n of Variable Depth, 

To evaluate the tides when the depth is a function of the latitude we must make 
use of the formulae of § 4. The method will be sufficiently illustrated by the 
computation of the solar semJ-diurnal tide for the case where 



or, where 



%/4a)%'^ = 2^0 + 3^- sin¥, 
h = (14,520 + 9680 sin^^^) feei 



denoting the co-latitude. 

Putting Ig/Aco^a^ = -^q, and making use of the numerical values found in § 15 for 

xl, y\, Ai, we obtain from (41) 



log ^l — 


3-()S907 


log ^^ = 


312901 


l0g7?l = 


2-17040 


log ^i = 


w. 375361 



log ^\ = n 3-06261 
log ^i — n 3-4369 

logTji =713-95767 
log 7?I = « 3-9962 



while, when -r-^., = -077, we obtain 



HI = + 0-020765 
iLj = - 0-039717 
m = - 0-052592. 






10 -~~ 



0*05787 
0-0606 



From these we deduce in succession, on neglecting M| 



IostIKTo = n 4-651 



'5 ^10 
log Ws 



r^ 4-2612 



logM 
logl^ 



log pll = 3-49214. 



But the first two of equations (40) give 






^^ I o __1Z^ 



72 



4-1632 
3-20891 






0„,Q, 



4wa- 



f 2^2 _ (1^2 _ |gJ2) q 






4 o 



4:(jd"ar 



m 



which, on solution, yield 



C! 



^-.9,^3 "t" ^ /<,.9^.o. ^ - 



Q2 



4 --^ 






4a)%' 



Q + 8 






m 



HI - mi 






(93 



©i 



4a)^a.^ ri lit' " M 



A; 



2 

'?3 



®| 



HI- mi 



7 + 8a?| 



kh 



®| 



4(u^a" 1tt(| 



mi 
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Oil substituting the numerical values for the quantities on the right, we obtain 

Q = - 2-9242 m. Q = 0-28546 C^| 

The remaining constants may now be computed from the formulae Cl/Gl = '^l/yjl, 
Cg/CI == Ml/y)t &c., and we finally obtain 

I = Qtll- 2-9242Po! + 0-28546P| - 0-00733P^ ~ 0-000147P| - 0-000007Pfo -.,..]. 

This makes the ratio of the height of the tide to that of the equilibrium tide at 

the equator 

- 3-6690. 

The tide will evidently be in the main inverted, the longest period of free oscilla- 
tion of the first class being in excess of twelve hours. 

As a further example, I have computed the series for C when the depth is given 
by the formula 

T¥ To ^il'^'^^i 



% _ _ 1 



that is, when the depth is 29,040 feet at the poles and shallows to 19,360 feet at the 
equator. This series is found to be 

C=mi- 2-3661P5 + 0'35649P| - 0'03953P^ + 0-00376P5 

- 0-000333Pio + 0'000029P'i. - 0-000002P?, + . . . / 

making the ratio of the tide to the equiUbrium tide at the equator 

- 3-4583, 

^, Ig 

If we put j~j^- = 3-Q, and replace X by 2cu in the period-equation, regarding this 

as an equation for /c, the largest root is found to be 

K = 21,765. 

Thus there will be a period of free oscillation coinciding exactly with twelve hours 
when 

h = (21,765 + 9680 sin^0 feet ; 

this formula makes the polar depth 21,765 feet, the equatorial depth 31,445 feet, and 

the mean depth 28,222 feet. 

let 
In like manner, when -^^ = ~ ifo, there will also be a period of free oscillation 
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agreeing exactly with twelve liours^ when the polar depth is oG/JTO feet^ th 
equatorial depth 26529O feet^ and the mean depth 29^517 feet. 



e 



§ 18. Diurnal Tides. 

It has been shovvni in § 14 that the diurnal tidal constituents whose periods are 
equal to a sidereal day will involve no rise and fall at the free surface when the 
depth of the ocean is uniform. This theorem will be rigorously applicable to the 
luni-solar diurnal constituent usually designated by the initial K^, while it may also 
be supposed to apply with a fair degree of accuracy to each of the solar diurnal 
constituents since the motion of the sun in his orbit is suiSficiently slow. There will 
however be an important lunar diurnal constituent for which the speed is 0'92700ca3 
in dealing with which we propose to take into account the difference between tlie^ 
period and a sidereal day. The method of computation is exactly similar to that 
used for the lunar semi-diurnal tides, and thus we find when hg/Aco^ci/' = -^q, 



^ = ©![- 0*07638Pi + 0-03543Pi - 0*00845P^ + 0'001207P^ 

- 0-000114P}o + O'OOOOOSPJo 

when hg/ioi^a^ =:= -^^o? 

I = ei[- 0169lP^i + 0'04738P1 - 0'00628P^ + 0'000480P^ 

— 0-000024P;o + O'OOOOOIPI^ 
when hg/ioxxt^ =^ ' 



J 



105 



'C-Qt},[^ 0'4145P.i + 0-06576Pi - 0'00464PJ + 0'000184P^ -- O'OOOOOoPlo + . . ^'\ 

and when hg/icD^a^^ = ^, 

^= ej[- l-4428P.^ + 0'123lPi - 0'00449PJ + 000009P^ «- 0'OOOOOlP}o + . . / 

It appears, then, that these tides will increase with the depth, and that they will 
be in the main inverted. For small depths the rise and fall will be small, but with 
a depth as great as 58,080 feet the tide will be in excess of the equilibrium»tide. 
The type will tend to approximate more and more closely to that represented by 
a second order harmonic alone as the depth increases^ 

So long as the depth is uniform the tidal constituents whose periods are rigorously 
ecpal to a sidereal day will never tend to become infinite, and consequently no 
period of free oscillation of the type of the diurnal tides can coincide exactly with 
one day. As o) diminishes however the shortest period of the second class, Avhich for 
the depths under consideration is longer than the period of the lunar diurnal 
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constituents, approximates to one day, and attains this as a limiting value when 
CO = 0. Hence, as co diminishes, or h increases, the largest root of the second class 
must pass through the value 0*92700, thus rendering one of the lunar diurnal 
constituents infinite. This accounts for the rapid increase in the coefficients in the 
series given above for these tides as h increases. 

The roots of the first class must however all be greater than unity, no matter 
how great the depth may be. Since they all decrease with the depth, they must 
approach finite limiting values greater than unity as the depth diminishes to zero. 

The tides of rigorously diurnal period will become infinite when the depth is 
variable if I assumes the value 



o 



SO that with this value of I we may anticipate that there will be a period of 
free oscillation exactly a sidereal day in duration. Now the above value of I will 
require that the surface of the solid earth should be rigorously spherical in order 
that the free surface of the ocean may be an equipotential surface under gravity and 
centrifugal force. It is easy to see why in this case a free oscillation of rigorously 
diurnal period must exist. For if the water be set in rotation as a solid body about 
an axis not rigorously coincident with the rotation-axis of the solid earth, and the 
form of the free surface be adjusted for equilibrium under centrifugal force about the 
new axis of rotation, there will be no forces acting which tend to modify this state 
of motion, and it will continue permanently, provided the system be free from 
friction. The motion of the water will be steady in space, but it will be oscillatory 
with a period of one day relatively to the solid earth. 

It is easy to verify that in the forced oscillations of rigorously diurnal period the 
motion of the water is ol* like character, involving no relative motion of the parts 
and being steady in space. In this case the axis about which the rotation of the 
water takes place lies in the plane containing the earth's polar axis and the dis- 
turbing body. 
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